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Experimental evidence for the sensitivity of the air-shower radio signal to the 

longitudinal shower development 
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We observe a correlation between the slope of radio lateral distributions, and the mean muon pseu- 
dorapidity of 59 individual cosmic-ray-air-shower events. The radio lateral distributions are mea- 
sured with LOPES, a digital radio interferometer co-located with the multi-detector-air-shower array 
KASCADE-Grande, which includes a muon-tracking detector. The result proves experimentally that 
radio measurements are sensitive to the longitudinal development of cosmic-ray air-showers. This 
is one of the main prerequisites for using radio arrays for ultra-high-energy particle physics and 
astrophysics. 



The investigation of ultra-high-energy cosmic rays 
works at the frontiers of particle physics as well as astro- 
physics. Only cosmic rays allow probing particle physics 
beyond the energy range of man-made accelerators. And 
only high statistics and accurate measurements of energy, 
direction and mass of cosmic rays will solve the mystery 
of their origin, i.e. how natural accelerators achieve en- 
ergies beyond 10^" eV (see references [ll, 0] for reviews). 
Since ultra-high-energy cosmic rays are too rare to be 
measured directly, they can only be measured indirectly 
by observing extensive air-showers of secondary particles. 
In particular, the longitudinal development of the air- 
showers is of interest for the two main physics goals: for 
particle physics, since the shower development depends 
on the interactions at ultra-high energies, and for astro- 
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physics, since the shower development also depends on 
the mass and type of the primary cosmic-ray particle. 

The longitudinal air-shower development can be stud- 
ied either with measurements of secondary air-shower 
particles, or by detecting electromagnetic radiation emit- 
ted by the air-shower or the traversed air. A method of 
the first category is the measurement of secondary muons 
with a muon-tracking detector (MTD), like it is per- 
formed at the KASGADE-Grande experiment [1,|1]: the 
directions of the muon tracks relative to the direction of 
the air-shower axis carry information on the longitudinal 
shower development (5|, |6[ . Gurrently, the most sophis- 
ticated methods for the longitudinal air-shower develop- 
ment fall into the second category. They are based on 
the detection of fluorescence or Gherenkov light emitted 
by air-showers, which both is limited to dark, clear and 
moonless nights - a disadvantage which can be overcome 
by measuring the radio emission of air-showers. Digital 
radio- antenna arrays can measure under almost any con- 
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FIG. 1. Shower cores of the 59 selected events measured si- 
multaneously with the east-west aligned LOPES antennas, 
the KASCADE-Grande particle-detector array and the muon- 
tracking detector (MTD). Due to a reconfiguration of LOPES, 
15 of the 30 antennas were only available for a part of the 
events. 



ditions and would increase the statistics by a factor of 
roughly 10 0- 

Until now, the sensitivity of radio measurements to 
the longitudinal air-shower development was predicted 
by simulations of the radio emission [g, but the exper- 
imental proof was still missing. We supply this proof 
by comparing the mean pseudorapidities of muons mea- 
sured by KASCADE-Grande with the slope of the radio 
lateral distribution measured with LOPES [3, M, Il3; a 
digital radio-antenna array co-located with KASCADE- 
Grande. LOPES features a precise time [llj| and abso- 
lute amplitude calibration [12|, and measures the radio 
signal between 40 and 80 MHz. Since the radio emis- 
sion is mainly of geomagnetic origin [3, \l^ [IJ i the radio 
signal is predominately east-west polarized. Thus, we re- 
strict this analysis to east-west aligned LOPES antennas, 
which for most events have a higher measurement quality 
than north-south aligned antennas. 

For the analysis we use 59 air-shower events with an 
estimated primary energy above 5 ■ 10^^ eV, which were 
measured between December 2005 and May 2009 with 
both the MTD of KASCADE-Grande as well as with 
LOPES (figure [1]). The MTD consists of 16 individual 
muon telescopes with three horizontal layers and one ver- 
tical layer of streamer-tube chambers, providing a total 
detection area of 128 m^ for vertical muons [J|. It tracks 
individual muons, and allows a reconstruction of the ra- 
dial angle p, and the tangential angle r to the shower 
axis [15| measured by KASCADE-Grande with an angu- 
lar resolution of about 0.3° each. We limit the analysis 



to events with a distance between the MTD and the air- 
shower axis from 160 m to 320 m. For most of the less 
distant events the MTD is hit by too many particles to 
reconstruct the individual muon tracks, and for most of 
the more distant events the radio signal is too weak. We 
apply several quality cuts on the KASCADE-Grande re- 
construction, e.g., a fiducial area cut to guarantee an 
accurate reconstruction of the shower geometry and pri- 
mary energy. Furthermore, we apply several cuts on the 
LOPES measurement accepting only events with a ra- 
dio signal clearly distinguishable from noise. Of those 
events passing these cuts we exclude 1 event recorded 
during a thunderstorm, and 19 events for which the fit 
of the exponential-lateral-distribution function fails (see 
reference [lO| for a description of the cuts). To elimi- 
nate scattered, low-energy muons from the analysis [5], 
we use only muon tracks with p < 8° and r < 2°. Our 
results are based on the henceforth used 59 events which 
pass all these cuts. The cuts naturally influence the pre- 
sented results on a quantitative level (e.g., the numbers 
resulting from the performed fits), but the qualitative re- 
sult, namely that there is a correlation between the mean 
muon pseudorapidity and the slope of the radio lateral 
distribution, does not depend on the specific values for 
the cuts. 

As observable for the longitudinal development of 
the air-showers we use the mean pseudorapidity of the 
muons, where the pseudorapidity of each muon is calcu- 
lated as —ln{y'T^ + p^/2) [15]. Due to the experimen- 
tal geometry, i.e. the confined muon-to-shower-axis dis- 
tance range, the larger the pseudorapidity of a muon from 
a given air-shower the larger is the production height 
along the shower axis. In principle, the present analysis 
could be performed directly with the muon-production 
heights. However, the uncertainty in reconstruction of 
the distance to the shower axis introduces an additional 
uncertainty to the production height. Thus, muon pseu- 
dorapidities can be measured with a higher precision, and 
the reported result has a higher significance for the pseu- 
dorapidities than for the production heights. For each 
shower, we average over the pseudorapidities of all muons 
(on average 23 ± 9 muon tracks per event). Measured at 
a certain distance to the shower axis, large mean pseu- 
dorapidities mean that the shower development started 
at a large distance to the detector. For more details, 
see references [la, [ij]. The uncertainties of the mean 
muon pseudorapidities have been estimated by propagat- 
ing the angular resolution of 0.3° of the individual muon 
tracks. For this we use the bootstrap method, a numer- 
ical Monte-Carlo method which takes into account the 
uncertainties of each individual data point [1^ . 

The radio lateral distribution is the maximum ampli- 
tude of the radio signal in the individual LOPES anten- 
nas as a function of the axis-distance of each antenna 
(figure [2]), where the air-shower radio pulse is identified 
by cross-correlation beamforming '19]. We determine the 
slope of the radio lateral distribution by fitting an expo- 
nential function with a slope parameter, where the fit 
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Exponential fit function: 

£{d) = ^jga exp(- a • (J - 100 m)) 




Result of fit: 

red. x' = 0.82 
amplitude at 100 m 8, 
15.6 ±1.9|jV/m/MHz 
slope parameter a: 

8.51 ± 1.35 km^ 
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FIG. 2. Example for a radio lateral distribution measured 
with LOPES. 



-D 
Q. 
ns 

O 

-a 

3 



2.5 




Result of line fit: 

function y = offset 



red. X 
offset 
slope 







slope 
1.37 
4.86 ±0.12 
-0.117±0.017km 



steeper 



10 



12 



14 



16 



LOPES lateral slope [1 /km] 



FIG. 3. Mean muon pseudorapidity and the slope parameter 
of the radio lateral distribution for the 59 events. 



gives us also the statistical uncertainty of the slope pa- 
rameter. From simulations and theoretical arguments 
Q we expect that the slope depends on the geometri- 
cal distance of the antenna array to the air-shower cas- 
cade, and thus is sensitive to the longitudinal develop- 
ment of the electromagnetic air-shower component and 
therefore to the primary mass. A steep slope would indi- 
cate that the shower developed at a small distance to the 
radio-detector array. Consequently, we expect an anti- 
correlation between the mean muon pseudorapidity and 
the radio-slope parameter, if this theoretical prediction 
holds true. This is what we indeed observe (figure [2|). 

The strength of the correlation can be estimated by 
calculating the correlation coefficient r which is (— )1 in 
case of a perfect (in the absence of measurement uncer- 
tainties), linear (anti-)correlation, and in case of no 
correlation. We calculate r to -0.64 ± 0.09 (fig.©. The 
uncertainty of r has been estimated with the bootstrap 
method and corresponds to a significance of 7. 1 tr for the 
correlation, i.e. r > is excluded with 7.1 cr. As a con- 
sistency check, the significance has also been estimated 



by fitting a line, minimizing an effective x^ taking into 
account x- and y-uncertainties of each data point [20j. 
Counting by how many sigmas the slope of the line is 
different from results in the same significance of 7.1 cr 
as for the bootstrap method. 

The found correlation is cither due to a causal link (the 
common sensitivity of both observables to the shower de- 
velopment), or due to common systematic effects and bi- 
ases of the measurement devices. In particular, we have 
checked that the correlation is not only a consequence of 
known detector effects related to the axis distance and 
the shower inclination (figure HJ. Moreover, possible ad- 
ditional effects due to the geomagnetic angle (= angle be- 
tween shower direction and geomagnetic field) have been 
investigated. Although the geomagnetic angle has a large 
influence on the radio amplitude, this is not true for the 
slope of the radio lateral distribution. Indeed, we observe 
no significant effect, neither for the slope parameter nor 
for the mean muon pseudorapidities. We have removed 
the systematic axis-distance effect from the correlation by 
correcting the measured mean pseudorapidity and slope 
parameter of the radio LDF according to a linear fit to 
the measured axis-distance dependence. After this cor- 
rection, the correlation is, as expected, slightly weaker 
(r — —0.62 ± 0.10), but still has a significance of 6.0 a. 

The systematic effect of the shower inclination can, 
in principle, be treated in a similar way, although the 
physics behind the effect is more complex. Even with- 
out any detector effects, we expect a correlation of both 
observables with the zenith angle, since the geometrical 
distance between the detector and the development of the 
air-shower cascade depends not only on the atmospheric 
depth, but also on the zenith angle 9. Without any model 
assumptions it is impossible to disentangle detector from 
physics effects. Thus, we correct both observables for the 
observed zenith angle dependence as we did for the axis- 
distance effect. As expected, the correlation (without 
axis-distance correction) gets weaker after the correction 
for the shower inclination (r ~ —0.56 ±0.11, significance 
of 5.3 cr). When iteratively correcting for both the axis- 
distance and zenith-angle dependencies, the result stabi- 
lizes after a few iterations, and there is still a significant 
correlation (figure [51 r = —0.44 ± 0.12, significance of 
3.7cr). 

We investigated whether \r\ < 1 indicates a non- 
linearity of the correlation between both observables. 
Therefore, we studied by how much \r\ would be reduced 
due to the measurement uncertainties of both observ- 
ables if the underlying relation were a perfect linear anti- 
correlation. In all cases (before and after corrections) 
the measured value of r is compatible with a linear re- 
lation between the mean muon pseudorapidities and the 
radio-lateral-slope parameter. Hence, more precise mea- 
surements or increased statistics will be necessary to test 
if the correlation is truly or only in first order linear. 

We conclude that the observed correlation which re- 
mains after the corrections, cannot be explained by 
known systematic detector effects. Consequently, the 
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Result of line fit: 

red. x' = 2.55 

offset = -0.76 ± 1.35 km 

slope = 31.9 ±5.7 km 2 
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FIG. 4. Dependencies of the two observables, mean muon pseudorapidity (left) and slope parameter of the radio lateral 
distribution (right), on the mean distance from the measured muons to the air-shower axis (top) and the shower inclination 
(bottom). The slopes of the fitted lines are used for correction of the observables. 
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FIG. 5. Correlation between the mean muon pseudorapidity 
and slope of the radio lateral distribution after correcting for 
the distance to the shower axis and the shower inclination. 



correlation is due to the air-shower development in the at- 
mosphere. The fact that the correlation becomes weaker, 
but is still present after zenith correction, is a strong 
indication that the slope of the radio lateral distribu- 
tion is primarily sensitive to the geometrical distance 



between the antenna array and the air-shower cascade, 
and not its atmospheric depth in g/cm^. The correlation 
is of particular interest, since the radio emission conies 
from the electromagnetic component of the air-shower - 
mainly due to the geomagnetic deflection of electrons and 
positrons [21| and the variation of the net charge dur- 
ing the shower development [22'| (see reference [23] for 
a recent comparison of radio-emission models). There- 
fore, the correlation is also experimental evidence for the 
strong link between the development of the electromag- 
netic and the muonic air-shower component, which itself 
is a signature of the hadronic component [13, [Ij, [25| . As 
a consequence, the radio lateral slope must be sensitive - 
as the mean pseudorapidity is - to all parameters which 
influence the shower development, in particular the en- 
ergy and mass of the primary particle, and the hadronic 
interactions. A discussion of those dependencies is, how- 
ever, beyond the scope of this initial analysis. 

To exploit the sensitivity of the air-shower radio emis- 
sion to the shower development for particle and cosmic- 
ray physics, quantitative parameters for the shower de- 
velopment like the atmospheric depth of the shower max- 
imum, Xmax, have to be reconstructed. This requires ad- 
ditional input by air-shower models or a cross-calibration 
with another technique, e.g., with fluorescence measure- 



ments of X^ax at the Pierre Auger Observatory [26| , or 
with air-Cherenkov measurements at the Tunka experi- 
ment (H- 

At LOPES, a radio reconstruction of Xmax has al- 
ready been performed using simulations as input 28], 
namely the simulated dependencies of Xmax on the ra- 
dio lateral slope [3l| and the radio wavefront [321 ■ The 
A^max precision of an individual LOPES event is in the 
order of lOOg/cm^, respectively 200g/cm^, when using 
the first, respectively the second method. The reason for 
this large uncertainty is the high level of anthropogenic 
radio noise at the LOPES site. Nevertheless, the simula- 
tions indicate that in a situation with negligible noise per- 
event precisions of better than 30 g/cm^ can be achieved, 
i.e. similar to the precision of the fluorescence and air- 
Cherenkov technique. To cross-check these simulation- 
based results, we have estimated the per-event Amax pre- 
cision from the correlation of figure [5l When converting 
the typical mean muon pseudorapidity of 4 at a zenith 
angle of 30° and a typical axis distance of 240 m to an 
atmospheric depth, the average per-event uncertainty of 
the radio lateral slope corresponds to a Amax uncertainty 
of about 115g/cm^, which thus is in the same order as 
the simulation-based results. 

Summarizing, a significant correlation has been found 
between the mean muon pseudorapidity and the slope 



of radio lateral distributions of individual air-showers. 
Since the sensitivity of the mean muon pseudorapidity 
to the longitudinal shower development is established 
[a,Hl3j the result serves as first experimental proof that 
also the radio lateral distribution is sensitive to the longi- 
tudinal air-shower development, which implies a sensitiv- 
ity to the primary mass of high-energy cosmic rays. Since 
the precision of LOPES is limited by human-made noise, 
dedicated experimental devices in regions with a low level 
of ambient noise are required to test whether the sensitiv- 
ity of radio measurements to the shower development can 
be used for a reasonable precise reconstruction of Amax 
and the mass of the primary cosmic-ray particles. 
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